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Cathode  reaction  models  for  solid  oxide  fuel  cells  with  proton  conducting  electrolyte  (H-SOFC)  are  pro¬ 
posed,  and  the  reacting  orders  for  each  elementary  step  with  respect  to  oxygen  and  water  vapor  partial 
pressure  are  calculated.  The  limiting  steps  of  cathode  reactions  are  investigated  with  Sm0.5Sr0.5CoO3_5 
(SSC)-BaCeo.8Sm0.203_5  (BCS)  composite  cathodes.  The  results  suggest  that  the  migration  of  protons  to 
TPBs  and  the  surface  diffusion  of  0~d  might  be  the  limiting  reactions  for  SSC-BCS  composite  cathodes  in 
wet  atmosphere,  while  the  oxygen  ions  transferring  into  electrolyte,  the  reducing  of  0 ac*  to  0~d,  and  sur¬ 
face  diffusion  of  0~d  might  be  the  limiting  reactions  for  SSC-BCS  composite  cathode  in  dry  atmosphere. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  with  proton  conducting  electrolytes  (H- 
SOFC)  have  attracted  much  attention  for  their  great  fuel  utilization 
efficiency  and  simple  fuel-recycling  instruments  [1-3].  In  H-SOFC, 
hydrogen  is  oxidized  at  the  anode  to  proton,  which  migrates  selec¬ 
tively  through  the  electrolyte  to  the  cathode,  and  undergoes  a 
half-cell  reaction  with  oxygen  to  produce  water.  The  overall  half 
reactions  at  the  anode  and  cathode  could  be  expressed  as  Eqs.  (1) 


and  (2),  respectively: 

H2  +  20*  -*  20H0#  +  2e' 

(i) 

40H o*  +  02  +  4e'  ->  2H20  +  40* 

(2) 

with  water  generated  at  cathode,  the  cathode  reactions  of  H-SOFC 
differ  a  lot  with  those  of  cells  with  oxygen  ion  conducting  elec¬ 
trolyte  (O-SOFC),  of  which  oxygen  is  reduced  by  the  cathode  to 
oxygen  ions  via  the  overall  half-cell  reaction: 

02  +  4e'  4-  2Vo##  ->•  20q  (3) 

The  difference  in  cathode  reactions  between  H-SOFC  and  0- 
SOFC  results  in  different  cathode  processes.  But  until  now,  few 
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works  have  been  reported  on  the  cathodic  processes,  especially  the 
limiting  steps. 

Uchida  et  al.  [4]  investigated  the  oxygen  partial  pressures  (Pq2) 
dependence  of  cathode  polarization  resistance  (Rp)  of  Pt  cath¬ 
ode,  and  proposed  the  elementary  reaction  steps  for  H-SOFC 
cathode.  Yamaura  et  al.  [5]  investigated  the  polarization  resis¬ 
tances  of  Lai_xSr*Fe03_(5  and  Pt  also  by  means  of  a  current 
interruption  method,  and  found  that  the  rate-determining  step  of 
La0.7Sr0.3FeO3_5  cathode  was  different  from  that  of  sputtered  plat¬ 
inum  cathode.  And  they  indicated  that  further  study  was  necessary 
to  specify  the  limiting  steps. 

In  this  work,  oxygen  reactions  at  Sm0.5Sr0.5CoO3_5 
(SSC)-BaCeo.8Smo.203_5  (BCS)  composite  cathodes  were  investi¬ 
gated  for  H-SOFC  with  BCS  electrolyte.  The  polarization  resistance 
of  SSC-BCS  electrode  was  investigated  by  A.C.  impedance  as 
function  of  oxygen  (Pq2)  and  water  vapor  partial  pressure  (Ph2o). 
respectively,  to  determine  the  rate-limiting  steps.  The  cathode 
reaction  mechanism  for  H-SOFC  was  proposed,  and  the  depen¬ 
dence  of  reacting  speed,  accordingly  the  polarization  resistance, 
on  Po2  and  Ph2o  were  calculated  for  each  elementary  step. 

2.  Experimental 

Symmetric  cells  with  Sm0.5Sr0.5CoO3_5  (SSC)-BaCe0.8Sm0.2O3_5 
(BCS)  composite  cathodes  were  fabricated  using  screen-printing 
technique  with  BaCe0.8Sm0.2O3_5  as  the  electrolyte.  The  proce- 
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Fig.  1.  Impedance  spectra  of  an  SSC-BCS  electrode  measured  at  500  °C  in  humidified 
N2-O2  atmospheres  with  various  oxygen  partial  pressures. 


dures  were  summarized  as  follows:  polycrystalline  powders  of  SSC 
and  BCS  were  prepared  by  a  glycine-nitrate  process  as  described 
in  details  previously  [3,6].  BCS  powders  were  pressed  at  220  MPa 
into  disks  about  13  mm  in  diameter  and  0.54  mm  in  thickness,  and 
then  sintered  at  1500 °C  for  5h  as  electrolyte  substrates.  SSC  and 
BCS  powders  were  mixed  at  a  weight  ratio  of  6:4  to  prepare  the 
print-ink  with  ethocel  and  abietyl  alcohol  as  organic  binder.  The 
composite  electrodes  about  40  p,m  in  thickness  were  formed  by 
screen-printing  the  ink  onto  each  side  of  the  electrolyte  substrate. 
The  samples  were  fired  at  1000  °C  for  5h  in  air  to  form  the  sym¬ 
metric  cells.  Pt  paste  was  applied  onto  the  electrode  surface  as  the 
current  collector. 

Symmetric  cells  were  tested  from  400  to  700  °C  with  a  home- 
developed-cell-testing  system  by  A.C.  impedance  method  using 
an  electrochemical  workstation  (IM6e,  Zahner)  under  wet  (~3% 
H20)  or  dry  02-N2  atmospheres.  Measurements  were  conducted 
with  the  oxygen  partial  pressure  ranging  from  0.05  to  1.00  atm  by 
mass-flow-meter  (D08-2D/ZM)  with  the  general  atmosphere  pres¬ 
sure  (dry  or  wet)  controlled  at  1.00  atm.  The  water  vapor  partial 
pressure  of  0.03, 0.06, 0.10, 0.15  and  0.20  atm  was  achieved  by  heat¬ 
ing  the  bubble  at  different  temperature.  The  frequency  range  for 
impedance  measurements  was  from  0.01  Hz  to  1  MHz.  Curve  fitting 
was  performed  using  the  Zview  software. 

3.  Results 

3.1.  Properties  of  SSC-BCS  composite  electrodes  under  wet 
atmosphere 

Fig.  1  presents  the  impedance  spectra  of  SSC-BCS  electrode 
measured  at  500  °C  in  wet  (~3%  H20)  02-N2  atmosphere  with 
oxygen  partial  pressure  ranging  from  0.05  x  105  to  0.97  x  105  Pa. 
Two  depressed  arcs  are  observed  in  each  spectrum,  a  small  high- 
frequency  arc  and  a  large  low-frequency  arc  with  relaxation 
frequencies  of  24  kHz  and  0.1  Hz  at  500  °C,  respectively,  implying 
that  there  are  two  rate-limiting  steps.  An  equivalent  circuit  com¬ 
posed  of  two  RQ  elements  (RhQh)  (PlQl)  are  proposed  to  resolve 
these  spectra,  where  R  represents  the  polarization  resistance,  Q 
represents  the  constant  phase  element,  and  the  subscripts  H  and 
L  correspond  to  the  high-  and  low-frequency  arc,  respectively.  At 
oxygen  partial  pressure  (Po2 )  of  0.97  x  105  Pa,  the  fitted  Rh  and  Ri  is 
2.07  and  5.23  £2  cm2,  respectively.  With  the  decrease  of  oxygen  par¬ 
tial  pressure,  PL  increase  largely,  while  RH  keeps  almost  constant. 
And  at  Po2  of  0.05  x  105  Pa,  Rh  and  Ri  are  2.17  and  10.78  Q  cm2, 
respectively,  suggesting  that  the  resistance  associated  to  the  low- 
frequency  arc  is  the  main  contribution  to  the  polarization  resistance 
at  500  °C.  The  values  of  capacitance  at  500  °C  are  about  0.5  pT  cm-1 
and  0.4  F  cm-1  for  the  high-  and  low-frequency  arc,  respectively. 

The  dependence  of  PH  and  PL  on  oxygen  partial  pressure  (P02 )  are 
presented  in  the  usual  form,  Rt  <  (Pq2  )-m,  where  m  is  also  deemed 
as  reaction  order  of  Pz  with  respect  to  Po2.  As  shown  in  Fig.  2,  PH 
is  almost  independent  on  oxygen  partial  pressure  over  the  entire 


Fig.  2.  Oxygen  partial  pressure  (Po2 )  dependence  of  high-frequency  resistance  (Rh, 
solid)  and  low-frequency  resistance  (Ri,  open)  measured  in  humidified  N2-O2  atmo¬ 
sphere  at  500  and  600  °C,  respectively. 

Po2  range  measured  with  fitted  mu  very  close  to  zero,  0.02  ±  0.01 
on  average,  suggesting  that  neither  atomic  oxygen  nor  molecular 
oxygen  is  involved  in  this  step.  While  for  PL,  a  linear  increase  of 
log(l/PL)  is  observed  with  the  increase  of  oxygen  partial  pressure, 
and  the  fitted  mL  is  close  to  1  /4  (0.24  ±  0.01 ). 

Arrhenius  plots  of  Rh  and  PL  at  different  oxygen  partial  pressures 
were  shown  in  Fig.  3  to  determine  the  activation  energy  for  react¬ 
ing  species.  Within  the  tested  temperature,  Rh  was  independent  on 
oxygen  partial  pressures,  and  the  simulated  activation  energy  was 
0.60  eV,  quite  close  to  that  of  proton  conducting  in  doped  BaCe03 
[7-10],  about  0.4-0.65  eV,  implying  that  the  reaction  correspond¬ 
ing  to  the  high-frequency  arc  might  be  proton  involved.  Although 
Ri  values  increased  with  the  decrease  of  oxygen  partial  pressures, 
its  activation  energies  are  about  1.18  eV  and  also  independent  of 
the  oxygen  partial  pressures.  The  high  activation  energy  of  Ri  sug¬ 
gests  that  Rl  decreases  more  rapidly  with  the  testing  temperature 
than  Rh.  At  temperature  below  600  °C  with  oxygen  partial  pres¬ 
sure  of  0.97  x  105  Pa,  Ri  exceeds  Rh  and  becomes  the  dominating 
contribution  of  electrode  resistance. 

The  dependence  of  Rh  and  PL  on  the  water  vapor  partial  pres¬ 
sure,  Ph2o>  are  also  investigated  with  Pq2  controlled  at  0.20  x  105  Pa, 
and  no  additional  arcs  besides  the  high-  and  the  low-frequency  arc 
are  observed  within  the  measured  water  vapor  pressure  range,  as 
shown  in  Fig.  4.  The  dependence  of  Rh  and  PL  on  Ph2o  are  quite 
different  with  those  on  Pq2.  At  500  °C  with  PH2o  increased  from 


Fig.  3.  Arrhenius  plots  of  Rh  (solid)  and  RL  (open)  measured  in  humidified  N2-O2 
atmosphere  with  oxygen  partial  pressure  Po2  =  0.20  atm  (■,  n),  0.5  atm  (a,  A),  and 
0.97  atm  (★,  ☆),  respectively. 
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Fig.  4.  Impedance  spectra  of  SSC-BCS  electrode  measured  at  500 °C  in  20%  02-N2 
atmosphere  with  various  vapor  partial  pressures. 


0.03  x  105  Pa  to  0.20  x  105  Pa,  PH  reduces  from  2.16  to  0.90  C2  cm2, 
while  Ri  keeps  almost  constant  at  7.89  C2  cm2,  suggesting  that  H20 
is  involved  in  the  steps  corresponding  to  the  high-frequency  arc  but 
not  in  that  to  the  low-frequency  arc. 

The  reaction  orders  with  respect  to  Ph2o,  n,  for  the  electrode 
reactions  are  determined  from  the  slope  of  R h  and  Ri  curves,  respec¬ 
tively,  as  shown  in  Fig.  5.  The  reaction  order  for  the  high-frequency 
reaction,  Hh,  is  close  to  1  \2  (0.46  ±  0.01 ),  while  the  reaction  order  for 
the  low-frequency  reaction,  nL,  appears  close  to  zero  (0.01  ±0.01). 
And  both  nH  and  nL  is  independent  of  temperature. 

Thus  the  dependence  of  high-frequency  resistance  and  low- 
frequency  resistance  on  Pq2  and  Ph2o  can  be  concluded  as  Eqs.  (4) 
and  (5),  and  their  detailed  steps  will  be  specified  by  later  discus¬ 
sion: 

Kh^S/h’o2  (4) 

Rl^Po2V4pS2o  (5) 

3.2.  Properties  of  SSC-BCS  composite  electrodes  under  dry 
atmosphere 

The  properties  of  SSC-BCS  composite  electrodes  are  also  inves¬ 
tigated  under  dry  02-N2  atmosphere,  in  which  BaCe0.8Sm0.2O3_5 
(BCS)  is  a  pure  oxygen  ion  conductor.  The  change  of  conduct¬ 
ing  species  from  proton  to  oxygen  ions  resulted  in  quite  different 
cathode  reactions.  Fig.  6  presents  the  impedance  spectra  for  the 
SSC-BCS  electrode  measured  at  500  °C  in  dry  02-N2  atmosphere 
with  various  oxygen  partial  pressures.  Unlike  the  impedance  spec¬ 
tra  tested  under  wet  atmosphere,  there  are  three  depressed  arcs 
in  each  spectrum,  noted  as  high-frequency  arc,  middle-frequency 
arc  and  low-frequency  arc,  respectively.  With  the  decrease  of  Po2 , 
the  length  of  middle-  and  low-frequency  arcs  increase  obviously, 


Fig.  5.  Vapor  partial  pressure  (Ph2o)  dependence  of  high-frequency  resistance  (fiH, 
solid)  and  low-frequency  resistance  (fiL,  open)  in  humid  20%  02-N2  atmosphere  at 
500  and  600  °C,  respectively. 


Fig.  6.  Impedance  spectra  of  SSC-BCS  electrode  measured  at  500  °C  in  dry  N2-02 
atmosphere  with  various  oxygen  partial  pressures. 


while  that  of  high-frequency  arc  is  the  same.  Detail  analysis  of  these 
spectra  was  conducted  using  the  equivalent  circuit  with  three  RQ 
elements  (R'H  Qh  ){R'm  Qm)(^l  Ql),  where  R'  presents  the  polariza¬ 
tion  resistance,  and  the  subscripts  H,  M  and  L  denote  the  process 
corresponding  to  high-,  middle-  and  low-frequency  arc,  respec¬ 
tively.  The  fitted  capacity  for  high-,  middle-  and  low-frequency 
arc  are  1.43  x  10-7,  1  x  10-4  and  0.704F  with  the  relaxation  fre¬ 
quency  130  K,  180  and  0.18  Hz,  respectively.  Tested  at  500  °C,  R'u , 
R'm  and  R'  were  1.14, 6.98,  and  5.26  C2  cm2,  respectively,  with  Pq2  of 
0.20  x  103  Pa. 

Fig.  7  presents  the  dependence  of  R'H ,  R'M  and  R'L  on  the  oxy¬ 
gen  partial  pressure  tested  at  500  °C  in  dry  N2-02  atmosphere.  As 
shown,  the  reaction  order  with  respect  to  oxygen  partial  pressure, 
is  close  to  0  (0.0  ±0.01),  3/8  (0.34  ±0.02)  and  1/4  (0.19  ±0.02)  for 
R'w  R'm  and  R'v  respectively.  The  activation  energies  were  also  inde¬ 
pendent  on  oxygen  partial  pressure,  about  1.28, 1.00  and  1.17  eV  for 
R'w  R'm  and  R'v  respectively.  The  activation  energies  for  R'u  is  almost 
twice  of  that  in  wet  atmosphere,  suggesting  the  different  reacting 
species  in  wet  and  dry  atmosphere.  It  should  be  noted  that  the  reac¬ 
tion  order  and  the  activation  energy  are  almost  the  same  in  both 
wet  and  dry  atmospheres  for  low-frequency  resistances,  suggesting 
similar  reacting  mechanism.  The  dependence  of  R'w  R'M  and  R'L  on 
Po2  can  be  summarized  as  Eqs.  (6)— (8): 


p° 

J02 

(6) 

,  D-3/8 
'  J02 

(7) 

d-1/4 

o2 

(8) 

logP0  /  Pa 

Fig.  7.  Oxygen  partial  pressure  (Po2 )  dependence  of  high-frequency  resistance  (i?'H), 
middle-frequency  resistance  (R'u),  and  low-frequency  resistance  (J^)  in  dry  N2-02 
atmosphere  measured  at  500  °C. 
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4.  Discussion 

4.1  Cathode  reaction  models  and  data  analysis  in  wet 
atmospheres 

The  cathode  reactions  of  H-SOFCs  are  more  complex  than  those 
of  O-SOFCs  due  to  the  formation  of  H20  at  the  cathodes.  And  the 
cathode  reactions  of  FI-SOFCs  might  consist  of  several  steps  includ¬ 
ing:  (1)  surface  dissociative  adsorption  and  diffusion  of  oxygen 
along  with  charge  transfer;  (2)  proton  migration  to  triple  phase 
boundaries  (TPBs);  and  (3)  the  formation  and  desorption  of  FI20. 
In  Uchita’s  [4]  reaction  models,  the  proposed  cathode  steps  mainly 
came  from  those  of  O-SOFC,  and  the  individual  steps  of  proton  react¬ 
ing  were  seldom  considered.  In  this  work,  detailed  reacting  steps 
for  the  cathode  of  H-SOFCs  are  proposed,  and  the  reacting  rate  for 
each  step  are  calculated  as  function  of  oxygen  water  vapor  partial 
pressures.  The  elementary  steps  are  expressed  as  Steps  1-8: 

Step  1 

02(g)  20ad 

Step  2 

°ad  +  e~  °ad 

Step  3 

°ad  °TPB 

Step  4 

0TPB  +  e  -*  0TPB 

Step  5 

OH* o  electrolyte  ^TPB  "F  0 o 

Step  6 

H^pb  +  0TPB  ->  0Htpb 

Step  7 

^TPB  ~F  0Htpb  — H20tpb 


r3  =  k3a0-d-k'3a0-m 

(ii) 

=  M0-pB  exp  (-=§)  -  k4a0^  exp  (2T) 

(12) 

r5  =  ^5aoH0electrolyte*  -  fc5aH+pB 

(13) 

r6  =  l<6aH+  a02-  -  k'6a0H- 

nTPB  UTPB  TPB 

(14) 

r7  =  /<7aH+a0H-pB  -  /4aH2oTPB 

(15) 

r8  =  k8  Qh2oTpb  ~  k'8Pu2o 

(16) 

where  r2-  is  the  general  reaction  rate  with  the  subscripts  denoting 
the  step  number;  /q  and  k!{  are  the  rate  constants  for  the  forward 
and  backward  reactions,  respectively;  aM  represent  the  activity  of 
appointed  species  marked  as  the  subscript;  and  R ,  T  and  F  is  the  gas 
constant,  absolute  temperature  and  the  Faraday  constant,  respec¬ 
tively.  The  symmetry  coefficient,  is  assumed  to  be  1/2.  E  is  the 
electrode  potential  generated  by  Nernst  equation,  as  shown  in  Eq. 
(17).  It  should  be  noted  that  the  protons  in  the  electrolyte  is  gen¬ 
erated  by  the  incorporation  of  water  into  the  BCS  electrolyte,  as 
shown  in  Eq.  (18),  and  thus  noH*0  electrolyte  can  calculated  by  Eq. 


(19): 

f  RT  \ 

Eeq  =  const  +  f  —  J  lnPo2 

(17) 

H20(g)  +  Oq  +  Vo"  20Ho* 

(18) 

/  1  \1/2pl/2 
“OH-o  electrolyte  ~  \kH  J  H2° 

(19) 

where  kH  is  the  rate  constant  for  the  balance  reactions  of  Eq.  (18). 

Thus  the  reacting  rate  of  Steps  1  -8  can  be  given  as  function  of 
p02  and  Pu2o-  For  example,  if  Step  7  was  assumed  to  be  the  rate¬ 
determining  step  and  Steps  1-6  and  8  were  in  virtual  equilibrium, 
the  rate  of  Step  7  could  be  calculated  by  figuring  out  aH+ ,  a0H-  and 
aH20TPB  • 

(1)  aH+:  Submitting  Eq.  (19)  into  Eq.  (13),  we  obtained  the  activity 
of  protons  at  TPBs: 


(2)  Qh2oTpb  :  From  Eq.  (16),  the  activity  of  adsorbed  H20  at  TPBs  can 
be  expressed  as  function  of  Ph2o>  as  shown  in  Eq.  (21 ): 


Step  8 

H20Tpb  H20(g) 


where  Steps  1-2  describe  the  dissociative  adsorption  of  oxygen  to 
O-;  Step  3  describes  the  diffusion  of  O-  species  along  the  SSC  sur¬ 
face  to  the  triple  surface  boundaries  (TPBs);  Step  5  describes  the 
transfer  of  protons  from  the  electrolyte  (OH*0  electrolyte ) t0  the  TPBs; 
Steps  6  and  7  describe  the  reaction  of  protons  and  oxygen  species; 
and  Step  8  is  desorption  of  H20  from  TPBs  into  gas  phase.  In  this 
model,  the  oxygen  dissociated  adsorption  and  diffusion,  Steps  1-4, 
are  similar  to  that  in  O-SOFCs  [11  ].  The  rate  equations  correspond¬ 
ing  to  Steps  1-8  can  be  calculated  by  r  =  rforward  -  rreVe rse,  and  the 
details  are  shown  as  follows: 


n  =  k,Po2  -  k\alad 


r2  =  k200ad 


-  k'2a0-  exp 

ad 


O) 

(10) 


aH2OrPB  =  PH20  (21) 

(3)  a0H-  :  The  activity  of  OH-  at  TPBs  can  be  calculated  from  the 
equilibrium  of  Eq.(9),  (10),  (11),  (12)  and  (14).  From  Eq.  (9),  the 
activity  of  dissociated  oxygen  (Oad)  can  be  obtained  as  Eq.  (22): 

“oa(i=(!|)  (22) 


Combining  Eqs.  (22)  and  (10),  we  obtain: 


Substituting  Eq.  (23)  into  Eq.  (11)  gives: 
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Substituting  Eq.  (24)  into  Eq.  (12)  gives: 


Substituting  Eqs.  (20)  and  (25)  into  Eq.  (14)  gives: 


Thus,  the  steady  state  reaction  rate  of  Step  7  can  be  calculated 
by  substituting  Eqs.  (20),  (21)  and  Eq.  (26)  into  Eq.  (15),  as  shown 
in  Eq.  (27): 


Thus  the  forward  and  the  reverse  reaction  rate  for  Step  7  can  be 
expressed  as  Eqs.  (28)  and  (29),  respectively. 

^forward  —  ^7^H20  (28) 


^reverse  -  l<7  Pu2° 


As  shown  in  Eq.  (28),  the  forward  reaction  rate  is  the  function  of 
Ph2o  and  independent  of  P0l .  Substituting  the  Nernst  equation,  Eq. 
(17)  into  Eq.  (29),  the  reverse  reaction  rate  is  obtained  as  Eq.  (30), 
which  is  also  the  function  of  Ph2o  and  independent  of  Pq2. 


Thus  the  reaction  order  for  Step  7  with  respect  to  oxygen  partial 
pressure  and  water  vapor  partial  pressure  are  0  and  1,  respectively. 

We  now  consider  Step  3  as  the  rate-determining  step  and  assum¬ 
ing  that  Steps  1, 2, 4-8  were  in  virtual  equilibrium.  The  reaction  rate 
of  Step  3  can  be  calculated  by  figuring  out  an-  and  an-  . 

ad  TPB 

(1)  a0-  :  Calculated  from  the  equilibrium  of  Eqs.  (9)  and  (10),  as 

ad 

shown  in  Eq.  (23): 

(2)  an~  :  Substituting  Eqs.  (20),  (21 )  into  Eq.  (15)  gives: 

TPB 

a°HT-PB = kV2  (i)  (kj)  (i)  p^°  (31) 


Submitting  Eqs.  (20)  and  (31),  into  Eq.  (14),  we  obtained: 


Table  1 

Elementary  cathode  reaction  steps  and  their  reaction  order  with  respect  to  oxygen 
partial  pressure  (m)  and  water  vapor  partial  pressure  (n),  for  solid  oxide  fuel  cells 
with  proton  conducting  electrolyte. 


Combining  Eqs.  (32)  and  (12),  we  obtained: 


We  obtained  the  reaction  rate  of  Step  7  by  substituting  Eqs.  (33) 
and  (23)  into  Eq.  (11): 


Substituting  the  Nernst  equation  into  Eq.  (34),  we  obtained: 


Thus,  the  dependence  of  r3  on  Po2  and  Ph2o  is  derived  as  Eq.  (36). 
And  the  reaction  rate  is  determined  as  1  / 4  and  0  with  respect  to  Po2 
and  Ph2o»  respectively. 

(36) 

And  the  reacting  speed  of  other  elementary  step  can  be  calcu¬ 
lated  in  the  same  way,  as  shown  in  Table  1,  where  m  and  n  refers  to 
the  reacting  order  for  each  step  with  respect  to  Po2  and  Ph2o>  respec¬ 
tively.  As  shown,  the  dependence  of  reacting  speeds  for  Steps  1  -4  on 
Po2  is  similar  to  those  in  O-SOFC  proposed  by  Kim  et  al.  [  11  ],  indicat¬ 
ing  similar  reaction  mechanism  for  oxygen  dissociated  adsorption 
and  diffusion.  It  also  should  be  noticed  that  Steps  5  and  6  have  the 
same  dependence  on  Pq2  and  Ph2o>  implying  that  they  cannot  be 
discerned  by  only  investigating  their  dependence  on  oxygen  and 
water  vapor  partial  pressure. 

4.2.  Data  analysis  for  SSC-BCS  in  wet  and  dry  atmosphere 

In  wet  atmosphere,  the  low-frequency  resistance  is  independent 
of  Ph2o»  while  proportionate  to  Pq2  with  the  reaction  order  of  1  / 4,  as 
shown  in  Eq.  (4).  As  referred  to  Table  1 ,  Step  3  has  the  same  depen¬ 
dence  on  Pq2  with  the  low-frequency  resistance,  suggesting  that 
the  low-frequency  arc  might  correspond  to  the  surface  diffusion  of 
0~d.  The  activation  energy  of  low-frequency  resistance  is  1.17  eV, 
close  to  that  in  SSC-LSGM  [12],  which  is  deemed  corresponding  to 
the  surface  dissociative  adsorption  and  diffusion  of  oxygen.  This  fur¬ 
ther  confirms  the  surface  diffusion  of  0~d  in  low-frequency  reaction. 
While  the  high-frequency  resistance  is  proportionate  to  Ph2o  with 
the  reaction  order  of  1/2,  and  is  independent  of  Pq2.  As  shown  in 
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Table  1,  both  Steps  5  and  6  might  be  the  corresponding  steps,  which 
describe  the  protons  transfer  from  the  bulk  of  electrolyte  to  TPBs 
and  the  reacting  of  protons  at  TPBs  with  O2-  to  form  OH-,  respec¬ 
tively.  Thus  the  specific  step  corresponding  to  the  high-frequency 
arc  cannot  be  discerned  simply  from  its  dependence  on  oxygen  and 
water  vapor  partial  pressure.  Yet,  considering  the  activation  energy 
of  high-frequency  resistance  of  0.6  eV,  Step  5  is  more  promising, 
since  the  generating  energy  of  OH-  is  usually  about  an  order  lower 
than  that  of  proton  migration  [13].  Thus  the  migration  of  protons  to 
TPBs  and  the  surface  diffusion  of  0~d  might  be  the  limiting  reactions 
for  SSC-BCS  composite  cathode  in  wet  atmosphere. 

In  dry  atmosphere,  the  low-frequency  resistances  have  almost 
the  same  the  reaction  order  and  the  activation  energy  with  that 
in  wet  atmospheres.  This  suggests  surface  diffusion  of  0~d  might 
also  be  the  corresponding  reactions  for  low-frequency  arc  in  dry 
atmosphere,  which  accords  with  that  for  O-SOFC  previously  [12]. 
Although  the  high-frequency  resistances  have  the  same  reaction 
orders  with  respect  to  Pq2  ,  their  activation  energy  is  1.28  eV,  about 
twice  of  that  in  wet  atmosphere,  suggesting  different  conducting 
species.  For  O-SOFC,  the  high-frequency  resistance  is  usually  related 
to  the  oxygen  ion  transfer  from  the  TPB  into  the  electrolyte  lattice, 
as  described  in  Eq.  (37),  because  of  its  independence  on  oxygen 
partial  pressure[ll,12].  This  differs  a  lot  from  that  in  H-SOFC,  sug¬ 
gesting  that  oxygen  ions  transition  might  the  corresponding  steps 
for  SSC-BCS  in  dry  atmosphere: 

OtFb+vO' '^0*  (37) 

While  the  reactions  order  for  middle-frequency  resistance  with 
respect  to  P02  is  3/8,  which  was  not  observed  in  the  spectra  in 
wet  atmosphere.  And  reducing  of  Oad  to  0~d,  as  shown  in  Table  1, 
is  proposed  to  the  corresponding  step  for  middle-frequency  arc. 
Thus  the  limiting  reactions  for  SSC-BCS  composite  cathode  in  dry 
atmosphere  should  be  the  oxygen  ions  transfer  from  TPBs  into  elec¬ 
trolyte,  the  reducing  of  Oad  to  0~d,  and  surface  diffusion  of  0~d. 

5.  Conclusion 

Reaction  models  for  solid  oxide  fuel  cells  with  proton  conducting 
electrolyte  (H-SOFC)  are  proposed  for  composite  cathodes  in  this 
work,  and  the  rate  of  each  elementary  step  is  calculated  as  function 
of  oxygen  and  water  vapor  partial  pressure,  as  shown  in  Table  1. 

Using  Sm0.5Sr0.5CoO3_5  (SSC)-BaCe0.8Sm0.2O3_s  (BCS)  compos¬ 
ite  cathode,  the  impedance  spectra  under  various  atmosphere  were 


measured  to  investigate  the  cathode  limiting  reactions  for  H-SOFC. 
Two  arcs  were  observed  for  SSC-BCS  under  wet  conditions.  The 
high-frequency  resistances  are  proportionate  to  (Ph2o) _1^2  and 
independent  of  Pq2  with  the  activation  energy  of  0.6  eV,  suggesting 
that  the  migration  of  protons  to  TPBs  might  be  the  correspond¬ 
ing  step.  While  the  low-frequency  resistances  are  proportionate  to 
(Po2  )-1/4  and  independent  of  Ph2o>  which  might  correspond  to  the 
surface  diffusion  of  0~d. 

While  in  dry  atmosphere,  where  BCS  is  a  pure  oxygen  ion 
conductor,  there  were  three  arcs  in  the  impedance  spectra.  The 
high-frequency  resistances  are  independent  of  Pq2  with  the  acti¬ 
vation  energy  of  1.28  eV,  which  might  correspond  to  the  oxygen 
ions  transferring  into  electrolyte.  The  middle-frequency  resistances 
are  proportionate  to  (Po2)-3^8,  and  thus,  the  reducing  of  Oad  to 
0~d  might  be  the  corresponding  step.  While  the  low-frequency 

resistance  is  proportionate  to  (Po2  )-1^4,  similar  to  that  in  wet  atmo¬ 
sphere,  which  might  correspond  to  surface  diffusion  of  0~d. 
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